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Abstract 

We explore the phase diagram of two-flavor QCD at imaginary values of baryon and isospin chemi- 
cal potentials, /ub and /^ijso, analyzing the thermodynamic potential of QCD analytically and that of the 
Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model numerically. QCD has no pion condensation 
at imaginary and ^iso, and therefore has discrete symmetries that are not present at real /j,b and /Uiso- 
The PNJL model possesses all the discrete symmetries. The PNJL model can reproduce qualitatively lattice 
QCD data presented very lately. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) as a fundamental theory on strong interaction is well de- 
fined, since it is renormalizable and parameter free. However, thermodynamics of QCD is not well 
understood because of its nonperturbative nature. In particular, QCD phase diagram is essential for 
understanding not only natural phenomena such as compact stars and the early universe but also 
laboratory experiments such as relativistic heavy-ion collisions. Quantitative calculations of the 
phase diagram from first-principle lattice QCD (LQCD) have the well known sign problem when 
the baryon chemical potential (/ib) is real; for example, see Ref. [1] and references therein. For 
later convenience, we use the quark-number chemical potential /iq = /iB/3 instead of fiB- So far, 
several approaches have been proposed to circumvent the difficulty; for example, the reweighting 
method [2], the Taylor expansion method [sl] and the analytic continuation from imaginary yUq to 
real fXq [|4|-l6|]. However, those are still far from perfection. 

As an approach complementary to first-principle LQCD, we can consider effective models such 
as the Nambu-Jona-Lasinio (NJL) model [ITl-llin and the Polyakov-loop extended Nambu-Jona- 



Lasinio (PNJL) model [112 



3311 . The NJL model describes the chiral symmetry breaking, but not 



the confinement mechanism. The PNJL model is constructed so as to treat both the Polyakov loop 



and the chiral symmetry breaking |113|]. 



In the NJL-type models, the input parameters are determined at /iq = and T > 0, where T 
is temperature. It is then highly nontrivial whether the models predict properly dynamics of QCD 
at finite /iq. This should be tested from QCD. Fortunately, this is possible at imaginary fiq, since 
LQCD has no sign problem there. 

Roberge and Weiss found [|34D that the thermodynamic potential i^QCD(^q) of QCD at imagi- 
nary chemical potential /iq = iT6q has a periodicity ^?QCD(^^q) = 'i^QCD(6'q + 27rA;/3), showing 
that J7QCD(^q + 2nk/3) is transformed into ^2QCD(^q) by the Z3 transformation with integer k. 
This means that QCD is invariant under a combination of the Z3 transformation and a parameter 
transformation 9q ^ 9q + 2k7r/3 [|30|l . 

q^Uq, A^^UA^-^ -t/g{dM)U-\ 6^ ^ 6^ + 27rk/3, (1) 

where U (x, r) are elements of SU(3) with U{x, 13 = 1 /T) = exp(— 2z7rA;/3)[/(a;, 0) and q is the 
quark field. We call this combination the extended Z3 transformation. Thus, f^QCD{dq) has the 
extended Z3 symmetry, and hence quantities invariant under the extended Z3 transformation have 
the RW periodicity [boll. 



At the present stage, the PNJL model is only a realistic effective model that possesses both the 
extended Z3 symmetry and chiral symmetry BOll . This property guarantees that the phase diagram 
evaluated by the PNJL model has the RW periodicity in the imaginary /Xq region, and therefore 
makes it possible to compare the PNJL result with LQCD data IJil^ quantitatively in the imagi- 
nary /iq region. Actually, the PNJL mode 
the vector-type four-quark interaction [Sl-. 



succeeds in reproducing the LQCD data by introducing 

n 



ion and the scalar-type eight-quark interaction [lion. The 



QCp phase diagram in the real /iq region is predicted by the PNJL model with the parameter 
set 113 in that reproduces the LQCD data at imaginary /Xq. The critical endpoint can survive, even if 
the vector-type four-quark interaction is taken into account. 



LQCD has no sign problem also at finite isospin chemical potential (/iiso) OSh . This is true 
for both real and imaginary isospin chemical potentials, as explicitly shown in Sec. HIl For later 
convenience, we use the "modified" isospin chemical potential /ij = yUiso/2 instead of /iiso itself. 



Very recently, LQCD data were measured at both real and imaginary /ii 1360 and also in the case 



case 112 



that both ut and /iq are imaginary [|370. The PNJL model has already been applied to the real /ii 



23|] . but not to the imaginary /ij case. 



In this paper, we explore the phase diagram of two-flavor QCD at pure imaginary values of /iq 
and /ii, by analyzing the partition function of QCD analytically and the thermodynamic potential of 
PNJL numerically. As the primary result, we will show that the pion condensation does not occur 
at imaginary fij and /iq and hence isospin and baryon number are conserved. As a consequence of 
this property, i^qco has higher discrete symmetries at imaginary /ii and /iq than at real /ii and Un- 
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The PNJL model possesses all the symmetries, and then the model reproduces LQCD data |13( 
qualitatively at imaginary /ii and /iq. Finally, the phase diagram at imaginary /ij and /iq is predicted 
by the PNJL model. 

In Secini it is shown at imaginary /iiso and /iq that no pion condensation takes place and then 
QCD has some discrete symmetries. A simple explanation of the PNJL model is made in Sec. Hill 
and numerical results of PNJL calculations are presented in Sec. UVl Section |V] is devoted to 
summary. 



II. DISCRETE SYMMETRIES OF QCD 



Roberge and Weiss showed the RW periodicity in the one-flavor case [34], assuming that 
baryon number is conserved. Extending their proof to the two-flavor case, we will prove that 
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^QCB{Oq, Oi) has some discrete symmetries at imaginary //q and //j. In this proof, we first as- 
sume that baryon number and isospin, i.e., w-quark and d-quark numbers, are conserved, but this 
assumption is confirmed to be true at the end of this section. 

The thermodynamic potential ^2QCD(^q! ^i) (per unit volume) is related to the partition function 
Z{9q, 9i) as ^?QCD = —T\n{Z)/V, where V represents the infinite volume we are thinking. The 
functional integral form of Z in Euclidean spacetime with time interval re(0,/3 = l/T)is 

Z = J DqDqDA cxp [-S] , 
S = Jd''x [g(7,D, - 74/i + mo)q + , (2) 

where q — {qu, qd)^ is the two-flavor quark field, mo = diag(m„, m^) is the current quark mass, 
and Dj, is the covariant derivative. We take the isospin symmetric limit of m„ = = rriQ. 

The chemical potential matrix fi is defined hy fi = diag(/i„, /x^) with the w-quark number 
chemical potential and the o?-quark one {jid). This is equivalent to introducing the baryon and 
isospin chemical potentials, jiB and jiiso, coupled respectively to the baryon charge B and to the 
isospin charge 73: 

/i = /^qTo + /^IT3 (3) 

with 

where tq is the unit matrix and Tj (i = 1, 2, 3) are the Pauli matrices in flavor space. Note that 

/ii is half the isospin chemical potential (iiiso). For later convenience, the dimensionless chemical 
potentials, 9y (y = u, d, q, I), are introduced by jiy = iTdy. 
Now, we transform the quark field q as 

q (exp[i6'„r//3]g„, exp[i6'dr//5]gd)'^ = exp[i6'qr//3] (cos [9it/I3]to + i sin[6'ir//3]r3)g. (5) 
This transformation leads Z to 

Z = / DqDqDAexp [-S] , 
S^Jd'x [qi^^D, + mo)q + \F'^,] (6) 

with the boundary conditions 

qu{x,/3) = -exp[i{9q + 9i)]qu{x,0), 

qd{x, ^) = - exp[i(^q - 9i)]qd{x, 0). (7) 
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Under the Z3 transformation, i.e., the first and second transformations of ©, Z keeps the same 
form as ([61) , but the boundary conditions are changed into 

qu{x, /3) = - exp[i(6lq + 9i - 27ik/3)]qu{x, 0), 

qdix, 13) = - exp[^(^^q -01- 2Txk/?>)]qd{x, 0). (8) 

The functional form of Q with the boundary conditions ([8]) means Z{6^ — 27rA;/3, ^i). Since the 
Z3 transformation corresponds to the redefinition of fields in the path integration, we can reach the 
equality 

z{d^,ei) = z{e^-2Tik/?>,ei). (9) 

Further, using ([7]), one can see that 

z{e^,e,) = z{9^,e, + 27r), (lo) 
z{e^ + 71,61) = z{e^,ei + 'K). (ii) 

In the isospin symmetric limit m„ = m^, Z is invariant under the interchange u d. This 
means that 

z{e^,e,) = z{e^,-e,). (i2) 

Furthermore, Z is invariant under charge conjugation, when 9^ and 9i are transformed as 9q — > —9q 
and 9i — )• —9i. This indicates that 

Z{9q,9i) = Z{-9q,-9i). (13) 

Equations (fT2l) and (fT3l) show that 

Z(0q,^^l) = Z(-^^q,^^l). (14) 

Thus, Z is 6'q-even and 6'i-even. The relations (fTTI) . (fT2l) and (fT4l) lead to new ones 

Z{9q, n±9,) = Z{9q + TT, ±^i) = Z(^q + 7r/3, ^i), (15) 

Z(7r/3 - ^q, 0l) = Z(^q - Vr/a, ^^l) = + TT, ^l) = 9l + n) . (16) 

The thermodynamic potential i^qcD = — Tln(Z)/V and the chiral condensate a = dilqcD / dmo 
have the same symmetries as Z in (|9l)-(fT6l). 
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Making the fermionic path integration in Q, one can get the determinant det A with A = 

'-)yDy — 74/i + rho- This determinant is real, since fi* = —fi and then [351 

(det A)* = det A^ = det(75Z\75) = det Z\ . (17) 

Further, A has an explicit form of 

det A = det [mil + {a ■ D - ^ij)\(y ■ D - fij)] 

X det [mil +{a-D- ^idl)\(r ■ D - /irf/)] , (18) 

where / is the 2x2 unit matrix and a-D = ID4 + ia-D. Each of the first and second determinants 
on the right-hand side of (fTSi ) is the square of a real number. Hence, det A is positive in the case 
(i) that both /iq and /ij are imaginary. 



Simi 



then [1351] 



arly, in the case (ii) that ji^ is imaginary and /ij is real, fi satisfies fi* = —TifiTi and 



(det Ay = det A^ = det(75riZ\ri75) = det Z\ . (19) 

This shows that det A is real. Furthermore, the determinant is given by 

det A = det [mp + {a ■ D - ^idl)\(y ■ D - fij)] 

X det [mil +{a ■ D- ^J)\(y ■ D - ^J)] . (20) 

This determinant is also the square of a real number and then positive. Thus, in both cases of (i) 
and (ii), LQCD has no sign problem. 

The Polyakov loop ^ and its Hermitian conjugate are defined as 

i> = —TtL, # = —TyL\ (21) 

with 

f/3 



r f 

L(x) = Pexp i / drA^i^jT) 
L Jo 



(22) 



where V is the path ordering and A4 = iA^. These are not invariant under the extended Z3 
transformation dH), so that their vacuum expectation values do not have the RW periodicity. We 
then introduce the modified Polyakov loop and its Hermitian conjugate, 

if = exp{ief)S, ij = exp{-ief)S^ (23) 
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for / = u,d. These are invariant under the transformation O. Their vacuum expectation values 
= i'^f) ^iid = i^f) have the same symmetries as Z in (|9l)- (fTT1) : note that is the complex 
conjugate of \Pf because Z is real. 

In the chiral limit, QCD has the chiral SUi^{2) x 5'f/R(2) symmetry when f^iso = 0. However, 
at /iiso 7^ this symmetry is reduced to Ui.j'i{l) x Ui.^ji(l), where I3 = is the third component 
of the isospin operator. Evidently, this symmetry can also be presented as f/13 (1) x [/AI3 (1), where 
Ui^{l) is the isospin subgroup and UAi^il) is the axial isospin subgroup. Quarks are transformed 
under these subgroups as g — )■ exp{iaT3)q and q — )■ exp{ia'y5Ts)q, respectively. In the case of 
= TTT'd > 0, only the Ui.^{l) symmetry survives. 

When QCD vacuum keeps the 11^(1) and Ui^{l) symmetries, the baryon charge B = V{q'j4q) 
is either zero or integer and the isospin charge I3 = ^(^74/3^') is also either zero or half-integer. 
In the partition function Z of the baryon- and the isospin-charge operator, q'j^q and q'jihq, 
appear through the form exp(2i 6^1 574/35 + iOc^q'^^q). Therefore, 9i and 9^ have periodicities (fTOl) 
and (fTT)) . Meanwhile, if the pion condensation occurs, the Ui^{l) symmetry is spontaneously bro- 
ken and hence the isospin charge is neither zero nor half-integer anymore. In this situation, QCD 
vacuum does not have periodicities (fTOl) and (fTT)) . We will then prove that the pion condensation 
does not take place at imaginary ii^^o- Son and Stephanov [35] show for real /ijso that the pion con- 
densation emerges when ||Uiso| > "^tt, where is the pion mass. For simplicity, we take /iq = 0, 
because the quark-number chemical potential does not break the f/13 (1) symmetry. Following their 
discussion in Ref. |35|] . we use the chiral perturbation theory that is applicable at yUjso smaller than 
the chiral scale (the p meson mass). The chiral Lagrangian for pion field S E SU(2) with finite 

/^iso IS 

£^jj = I^TtV^EVuS^ I^ReTrr (24) 

with flavor covariant derivatives 

VoS = - ^ {tsE - rrs) , 
VoS^ = doS^ + ^ (rtrs - rg^t) , (25) 
where is the pion decay constant. In the effective theory, the condensate IJ is described by 

IJ = tq cos a + iTi sin a. (26) 
The tilt angle a is determined by minimizing the vacuum energy (the static part of £cfr) 

C% = [{x -af-l- a'] (27) 



with X = cos a and a = (m^//iiso)^. Here, the static part has been obtained by inserting (l26l) 
into (|24l) . Noting that — 1 < s < 1, one can find for real yUiso that the static Lagrangian becomes 
minimum at x = 1 (a = 0) when a > 1 (niso < m^) and at a; = a (a = arccos(m7r//Uiso)^) when 
a < 1 (piso > ^n) 1351]. The fact that x = 1 and then = tq at yUjso < rriir means that the pion 
condensation does not take place there. 

As expected from (fT2l) . the static Lagrangian is /iiso-even and then a function of yU^^. Hence, 
the static Lagrangian with imaginary isospin chemical potential /Xiso = ii^ is given by substituting 
iu for /iiso in (ITT]) : 

Cf, = -^-^[ix + br-l-b^] (28) 

with b = (m,r/ i^Y- This static Lagrangian is minimum at x = 1 for any value of u. Therefore, the 
pion condensation does not occur at imaginary fiiso- The PNJL model can reproduce this property, 
as shown in Sec.HlIl 

The absence of the pion condensation at imaginary niso can be understood intuitively as follows. 
For real fiiso, the Bose-Einstein distribution function has an infrared divergence at /iiso > "^tt- This 
induces the Bose-Einstein Condensation, that is, the pion condensation. For imaginary /ijso, such 
a divergence never happens and then no pion condensation occurs. 

Putting X = 1 in (l28l) . one can obtain 

C% = -if.^nf. (29) 

Thus, in the limit T — )■ 0, the static potential (the thermodynamic potential) is independent of 
imaginary /iiso. The PNJL model can reproduce this property, as shown later. 



Ill, PNJL MODEL 

The two-flavor PNJL Lagrangian in Euclidean spacetime is 

£ = g(7,D, - 74/i + mo)g + G's[(gg)' + (g^Ts^g)'] - U{$[A],$[AY ,T), (30) 

where D„ = d„ - iA^. The field A" is defined as A" = gAl^S^^ with the gauge field A^, 
the Gell-Mann matrix and the gauge coupling g. In the NJL sector, Gg denotes the coupling 
constant of the scalar-type four-quark interaction. The Polyakov potential hi, defined in (l38l) . is a 
function of the Polyakov loop <P and its complex conjugate <P*. In the case of nio = ni = 0, the 

8 



PNJL Lagrangian has the SUl{2) x SUk{2) x U^{1) x SUc{3) symmetry. In the case of mo 7^ 
and fii ^ 0, it is reduced to U^{1) x U^{1) x 5*^4(3). 



In the Polyakov gauge, L can be written in a diagonal form in color space U3 



I = gi/3(<A3A3+08A8) ^ diag(e^'^'^^e*'^'^^e^^'^^), 



(31) 



where (l)a = h + </'8/V3, h = -h + </>8/y3 and </)c = -((/)a + = -2</)8/y3. The Polyakov 
loop ^ is an exact order parameter of the spontaneous Z3 symmetry breaking in the pure gauge 
theory. Although the Z3 symmetry is not exact in the system with dynamical quarks, it still seems 
to be a good indicator of the deconfinement phase transition. Therefore, we use ^ to define the 
deconfinement phase transition. 

The spontaneous breakings of the chiral and the UiSV) symmetry are described by the chiral 
condensate a = (qq) and the charged pion condensate [|22n 



(32) 



where t± = (ri ± iT2)/V2. Since the phase ip represents the direction of the Ui.^{l) symmetry 
breaking, we take Lp = for convenience. The pion condensate is then expressed by 



TT = {qi-f^nq). 



(33) 



The mean field (MF) Lagrangian is obtained by i22|] 



-^^MF = qiluD,. - 74/i + Mtq + Ni'j5Ti)q 



(34) 



where M = mo — 26*80" and = — 2Gs7r. Performing the path integral in the PNJL partition 
function 



ZpNJL = / DqDq exp 



- / d^xC 



MF 



(35) 



one can obtain the thermodynamic potential i7 (per unit volume), 



n = -rin(ZpNjL)/v^ 

i=± 

1 



(2vr)= 



3i5^(p) + ^ in [1 + 3(^ + (P*e-^^-"(P))e-'^^«"(P) + e'^^^^'^P)] 



+ ^ in [1 + 3(^* + <l>e-P''tiv)y-PEtiv) + ^-WEti.p)^^ ^q^'' + 7r^]+U (36) 



with 



(37) 



-^±(p) = -^±(p) ± /^q and -E'(p) = a/p^ + M^. Obviously, Q does not have discrete symmetries 
([To]) and when tt ^ 0. 

On the right-hand side of (|36l) . only the first term diverges, and it is then regularized by the 
three-dimensional momentum cutoff A 11131 llTH . The parameter set, A = 631.5 MeV, = 
5.498 [GeV^^] and mo = 5.5 MeV, can reproduce the pion decay constant /^r = 93.3 MeV and 
the pion mass = 138 MeV at T = [10]. We then adopt these values for A, Go and mo. We 



use U of Ref. il8|] that is fitted to LQCD data in the pure gauge theory at finite T [|38 



-'^^<P*<P + h{T) ln(l - + 4(^3 + ^*^) - ?>{<P<P*f\ 



(38) 
(39) 



where parameters are summarized in Table I. The Polyakov potential yields a first-order decon- 
finement phase transition at T = Tq in the pure gauge theory. The original value of Tq is 270 MeV 
determined from the pure gauge LQCD data, but the PNJL model with this value of Tq yields 
somewhat larger value of the pseudocritical temperature at zero chemical potential than the full 



LQCD simulation 



y, 



4l|l predicts. Therefore, we rescale Tq to 212 MeV [|3l|l . 



ao 


ai 


02 


^3 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential sector determined in Ref. ilSQ. All 
parameters are dimensionless. 



The classical variables X = <1>,<^*, a and vr satisfy the stationary conditions, 

dn/dX = 0. (40) 

The solutions of the stationary conditions do not give the global minimum of H necessarily. There 
is a possibility that they yield a local minimum or even a maximum. We then have checked that 
the solutions yield the global minimum when the solutions X(9q, 9i) are inserted into (l36l) . 

Now we numerically confirm that the pion condensation does not occur at imaginary fij. For 
simplicity, we set /Xq = 0, since the quark-number chemical potential does not break the Ui^{l) 
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symmetry. For this purpose, we search for the potential minimum by varying <P* and a with tt 
fixed. The potential surface H^ir) thus obtained is a function of 7f and drawn in Fig.[Tl where T is 
taken to be 175 MeV. Three cases of = 0, 7r/2 and tt are represented by the solid, dashed and 
dotted curves, respectively. For the three cases, the global minimum is always located at vf = 0. 
The curvature around the minimum becomes large as 9i increases. This means that the vacuum 
becomes more stable for larger 9i. 

-18 
-20 

^ -22 

IG 

-24 
-26 

-0.05 0.05 

71 



1 




6|=0 '- 






nil 






ji 








1 1 1 



Fig. 1: Potential surface as a function of vr at T = 175 MeV and = 0. The solid, dashed and dotted 
curves denote three cases of Q\ = 0, 7r/2 and vr, respectively. 

Therefore, we can set vf = 0. In this situation, the transformation ([5]) reduces £mf of (|34|) to 

Cmf = q{luD, + Mro)g - G.a'' - U (41) 

with the boundary conditions (|7]). Note that this procedure breaks down if vf ^ 0, since the 
operator qi'j^Tiq is not invariant under the transformation ([5]). Following Sec. UIl one can show 
that the thermodynamic potential f] has the same symmetries as Z in (|9l)-(fT6l). This statement is 
proven below more explicitly. 

Under the fact that vf = 0, i7 is reduced to a simpler form 



f=u,d • 



(27r)n 13 ' ^ 

+G,a^ + U. (42) 



where £'^(p) = E{i))±fif = E{p) ±i9f / f3. Obviously, 12 has discrete symmetries (flOl) and (fTTI). 
In the limit of T = 0, on the right-hand side of (l42l) the first term including 3£'(p) and the term 
Gscr"^ + U survive, and hence f] has no /iq and fii dependences there. 
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The thermodynamic potential i7 of (l42l) is not invariant under the Z3 transformation, 

although W of (1381 ) is invariant. Instead of the Z3 symmetry, however, J7 is invariant under the 
extended Z3 transformation, 

g±ie, ^ g±ie,g±i2|t^ ^ ^ ^e-^^, ^* ^ ^*e*^. (44) 

This is easily understood as follows. It is convenient to introduce the modified Polyakov loop 
\Pf = e^^f<P and ^I'j = e~*^^<P* that are invariant under the transformation (|44|) and have the same 
symmetries as Z in (fTOl) -(fTTI). The extended Z3 transformation is then rewritten into 

^±ie^ ^ ^±ie^^±i^^ ^^^^^^ ^}^^}, (45) 

and i? is also into 

= -2 ^ j ^[3i5;(p) + ^ln [1 + 3!^ye-^^(P) + 3<^;e-2'^^(P)e3^^^ + e-^^'^^^P^e^^*' 

+ - In [1 + 3!Z>;e-^^(P) + 3>Z>.e-2/^^(P)e-3^^^ + e-3/3^(p)e-3^^J 
/3 

Obviously, Q is invariant under the extended Z3 transformation (l45l) . since it is a function of 
only extended Z3 invariant quantities, e^*^^ = e^^^^e^^*^' (+ for u-quark and — for d-quark) and 
X(= a). The explicit dependence appears only through the factor e^'^^'^ in (|46l) . Hence, 

the stationary conditions (|40l) show that X = X(e'^*^'i). Inserting the solutions back to (|46l) . one 
can see that i7 = i7(e'^*^''). Thus, X and i? have the RW periodicity, 

9rrk 97rk 

m + ^) = m)> f2{e^ + ^) = ^{e^), (47) 



-GsfT^ + U. (46) 



and then 



^(^q + ^) = e-*2W3^(^j_ (48) 



The thermodynamic potential i7 of (1461) is invariant under the transformation 9i — > —9i, indi- 
cating that n is 6'i-even. The thermodynamic potential is also invariant under the — t- —9q 
transformation, if \l/f is replaced by This means that the solutions of the stationary condition 
(|40|) satisfy 

^;(0q) = ^;(-^q), (49) 
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indicating that i? is ^q-even. Furthermore, f2 of (l46l) satisfies the symmetries (flOl) and (fTTI) . These 
properties, together with the RW periodicity, guarantee that H of PNJL has the same symmetries 
as Z of QCD in (|9l)-(fT6l). The symmetries (l9l)-(fT6l) are visualized by numerical calculations in Sec. 
HYl 

Particularly at 9i = tx 12, Q has aperiodicity of 7r/3 in 6'q, because taking Qi to 7r/2 in (15) leads 

to 

Q{e^, 7r/2) = Q{e^ + 7r/3, 7r/2). (50) 

As shown in (|42|) . 12 is a sum of the thermodynamic potential Qui^u) for w-quark and that Qdi^d) 
for d-quark, i.e., Q = f2u{Qu) + ^didd), and f2f{0f) is a periodic even function of 3^^/. Hence, Qf 
can be expanded by cos(3A;6'/) with integer k. We then have 

i7 = ^ flfc [cos(3A;^„) + cos(3A;^d)] . (51) 

At lower temperature such as T < 2Tc, where Tc is the pseudocritical temperature of the decon- 
finement transition at /Xq = IM = 0, the coefficients [ak] of the expansion have the property that 



the flfc with A; > 2 are small [1320. In particular when Oi = n /2, Q is reduced to 



Q^2aQ + ai [cof^{W^ + 37r/2) + cos(3eq - 37r/2)] = 2ao (52) 

for any ^^q. Accordingly, when Oi = n/2, Q has a periodicity of 7r/3 in ^^q, but the dependence is 
quite weak. This property is also visualized by numerical calculations in Sec. |IVl 

IV, NUMERICAL RESULTS 
A. 0q dependence 

^q dependence of Q, the quark number density = —dn/d(iT9q) and the isospin number 
density ni = —dQ / d{iT6i) is investigated in this subsection. The thermodynamic potential i? is 
real and ^^q-even, so that riq and rii are pure imaginary, riq is ^^q-odd and ^^i-even. rii is ^^q-even and 
^i-odd. 



As for 6i = 0, it is known that, at temperature above Trw = = 190 MeV MM, dQ/dO^ 

is discontinuous at 6'q = 7r/3 mod 27r/3; note that Tc = 173 MeV in the present PNJL calculation. 
This discontinuity is called the RW phase transition. At such higher temperatures, three Z3 vacua 
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emerge alternatively in variation of 9^, that is, the first vacuum appears in the region (I) — vr/S < 
9q < 7r/3, the second one in the region (II) n/3 < 9q < n and the third one in the region (III) 
— TT < < — 7r/3. As a result of this mechanism, dQ/d9^ becomes discontinuous at boundaries 
of the three regions 11331 13411 . The charge conjugation is an exact symmetry on the boundaries. It 
is preserved below T-^, but spontaneously broken above Trw MM- 
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Fig. 2: dependence of (a) j7/r^, (b) lm[nq]/T^ and (c) lm[ni\/T^ at T = 175 MeV. Three cases of 
= 0, 7r/2 and vr are represented by solid, dashed and dot-dashed curves, respectively. 

Now, we consider T = 175 MeV as a typical temperature below Tr;vv- Figure [21 presents 9^ 
dependence of il/T^, the imaginary parts Im[nq/T'^] and Im[ni/T^] for three cases of 9i = 0, 7r/2 
and 71. These quantities have the RW periodicity and are smooth at any ^^q, as expected. Further, i? 
and Til are ^^q-even, while Uq is ^^q-odd. In the case of = 7r/2, J7 is almost constant and Im[?2q] 
is then nearly zero, as expected from (|52|) : precisely, they have a periodicity of n/3, but the 
dependence is quite weak. Meanwhile, Im[rii] is zero when = and tt, because it is ^^j-odd and 
satisfies (fTTI) . As for the case of = 7r/2, Im[r?,i] has the RW periodicity clearly. 

Figure [3] shows the same quantities as Fig.[2l but its temperature is T = 250MeV higher than 
Trw The RW periodicity is seen also in this figure. In the case of 9i = 0, il and ni have cusps 
at = 7r/3 mod 27r/3, while is discontinuous there. This discontinuity means the RW phase 
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Fig. 3: 6lq dependence of (a) f?/T^, (b) lm[nq\/T^ and (c) lm[ni]/T^ at T = 250 MeV. Three cases 
of = 0,7r/2 and vr are taken. In the panel (c), the sohd and dot-dashed lines agree with the x axis. 
Definitions of curves are the same as in Fig.[2l In the insets, these quantities at = 7r/2 are magnified. 

transition. When 6i = 7c/2, is almost constant, as expected from (|52|) . and Im[nq] is tiny 
everywhere. In the insets where H and Im[nq] at 9i = n/2 are magnified, as expected from (fT6l) . 
f2 and Im[ni] have cusps at 6'q = mod 7r/3, while Im[?T,q] is discontinuous there. As for the case 
of 6*1 = 7r/2, thus, the RW phase transition occurs at = mod n/3. Equation (fT6l) yields a 
relation 



As a consequence of this symmetry, in Figs. [2] and [3l the dot-dashed curves are obtained by shifting 
the corresponding solid curves by vr/S in the 9q direction. 

The discontinuity between the right- and left-hand limits of Im[nq(^q)] as 9q approaches vr/S, 
i.e., Im[nq(+7r/3) — riq(— 7r/3)], decreases as 9i increases from and disappears at 9i = n /2 + 
5{T), as shown later in Fig. IT^b) and[8te). Here, 5{T) numerically obtained is a small number 
depending on T weakly: 



f]{9q-n/3,9i + n) = Q{9q,9{). 



(53) 



5{T) = 0.00016 X (T - 250) 



(54) 
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for T > 212 MeV. Since the discontinuity of Im[nq(^q)] means the RW phase transition, 9i = 
7r/2 + 5{T) represents a location of an endpoint of the RW phase transition. Further discussion on 
the endpoint is made later in Sec. IIVEI 

For simplicity, our discussion begins with the case of T = 250 MeV, since 5(T) = there. 
Figure |4] (a) presents dependence of il/T'^ for five cases of = 0, vr/S, 7r/4, Stt/S and it/ 2. 
These results show that the RW phase transition occurs at ^^q = vr/S mod 2tt/3 when < 9i < 7r/2. 
Figure |4](b) represents the location of the RW phase transition in 9q-9i plane by solid lines. As 
mentioned above, when < 9i < n/2, the RW phase transition occurs at ^q = n/3 mod 2tt/3. 
This RW phase transition is also seen at —n/2 < 9i < 0, because H is ^^i-even. Furthermore, (|53l) 
indicates that the RW phase transition occurs also at 6'q = mod 27r/3 when n/2 < 9i < 37r/2. 
This is clearly seen in Fig.lH 



(b) 



01234 01234 

eq/(7i:/3) eq/(7i:/3) 
Fig. 4: (a) 6'q dependence of fl/T'^ for five cases of 6*1 = (solid curve), vr/S (thick dashed curve), 7r/4 

(thin dashed curve), Svr/S (dotted curve) and 7r/2 (dot-dashed curve), (b) Phase diagram in ^q-^i plane. The 

solid lines represent the RW phase transition. For both the panels, T = 250 MeV. 

For other T larger than 212 MeV, 5{T) is not zero. This makes the situation a bit more com- 
plicated. Following the logic mentioned above, we can find that the RW phase transition occurs at 
9q = n/3 mod 2n/3 when -n/2 - 6{T) <9i <n/2 + 6{T) and also at = mod 2n/3 when 
n/2 - 6(T) < 9i < 3n/2 + 6(T). This behavior in the vicinity of 9i = n/2 is confirmed in Fig. [5] 
that presents the phase diagram in 9q-9i plane at (a) T = 220 MeV and (b) T = 300 MeV. Note 
that 5{T) is negative in panel (a), but positive in panel (b). 
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(b) 



" 1 2 3 4 " 12 3 4 

ej(7i/3) ej(7r/3) 

Fig. 5: Phase diagram in 9q-9i plane at (a) T = 220 MeV and (b) T = 300 MeV in the vicinity of Oi = tt/2. 
The solid lines represent the RW phase transition. 

B. 9i dependence 

01 dependence of (2, and rii is investigated in this subsection. Equations (fTOl) and (fT2l) lead 
to a relation 

f2{e^, n-ei) = r?(^q, - vr) = f2(^q, + vr). (55) 

Thus, 01 dependence of i? is symmetric with respect to the axis 9i = vr. Differentiating (155]) with 
respect to ^^q, one can see that ^q-odd quantities such as have the same symmetry as 6'q-even 
ones such as (2: 

nqiOq, 71-61) = n^iOq, 9i - tt) = n^{9^, 9i + -k). (56) 

In contrast, differentiating (1551 ) with respect to 0q leads to the fact that the Oi dependence of the 
6'i-odd quantities such as rii is asymmetric with respect to the axis 61 = tt: 

-ni(^q,7r-0i) =ni(^q,0i + 7r). (57) 

Taking ^^q to 7r/6 in (fT6l) . one can find 

n{n/6,9i) = n{n/6,9i + n), ni{7i/6,9i) = ni{n/6,9i + tt), (58) 

indicating that ^^q-even quantities such as and nj have a periodicity of tt in 9i when 9^ = tt/6. 
Similarly, differentiating (fT6l ) with respect to 9q and setting 9^ to n/Q, we can get 

nq(7r/6, ^i) = -nq(7r/6, 9i + n). (59) 
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Fig. 6: 01 dependence of (a) Q/T"^, (b) lm[nq\/T^ and (c) lm[ni]/T^ at T = 175 MeV. Three cases of 
0q = 0, 7r/6 and 7r/3 are taken. In panel (b), the solid and dot-dashed lines agree with the x axis. 

Thus, ?2q(7r/6, 6'i) has an anti-periodicity of n in 9i, that is, the sign of Uq is changed by the 
transformation Oj ^ Oj + n. These properties of (|55])-(l59l) are seen below in Figs. [6]and|7l 

Figure [6] presents 9i dependence of Q/T'^, Im[nq]/T'^ and Im[ni]/T'^ at = 0,7r/6 and 7r/3 
for the case of T = 175MeV that is just above Tc = 173MeV and below Trw The quantities i? 
and Im[r?,q] are symmetric with respect to the axis 9i = vr, while Im[r?,i] is asymmetric with respect 
to the axis, as predicted by (l55l) - (l57l) . These are smooth everywhere in 9i, and have a periodicity 
of 2n for all ^^q. For 9q = n /6, Im[ni] (Im[nq]) has a periodicity (anti-periodicity) of tt in 9i, 
as expected from (|58] ) and (|59l ). Below T^w, Im[nq] is smooth at any ^^q. Hence, ^^q-odd quantities 
like Im[raq] are zero at ^^q = and n/3 mod 2n/3. In panels (a) and (b), as a result of the property 
of (|52l) . all curves almost meet at 9i = 7r/2 and 3n/2. As predicted by (|53l) . in all panels, the 
dot-dashed curve for the case of ^^q = 7r/3 is obtained by shifting the solid one for the case of 
6'q = by TT in the 9i direction. 

Figure |7] shows the same quantities as Fig. [6l but T is taken to be 250 MeV as an example 
of temperature above Trw Again, i? and Im[nq] are symmetric with respect to the axis 9i = tc. 



18 




0.5 1 1.5 2 0.5 1 1.5 



2 -(c) \/e 

7t/3 




0.5 1 1.5 2 

Fig. 7: 6*1 dependence of (a) Q/T"^, (b) Im[nq]/r3 and (c) Im[ni]/T^ at T = 250 MeV. Three cases of 
9q = —0, tt/6 and 7r/3 — are taken. 

while Im[ni] is asymmetric with respect to the axis. All the quantities have a periodicity of 2% 
in 9i for all 9^. For 6'q = n/6, Q and Im[ni] have a periodicity of tt in 6i, while Im[?7,q] has an 
anti-periodicity of vr in ^^i. Hereafter, we consider the right-hand (left-hind) limit of f{x) as x 
approaches a and denote it by f{x) |a;=a±o- As predicted by (l53l) . the dot-dashed curve for the case 
of = 7r/3 — is obtained by shifting the solid one for the case of 6'q = — by tt in the 6i 
direction. All curves almost meet oX. Oi = n/2 and Sir/ 2. 9i-even quantities such as H and riq 
have cusps at 6'q = 7r/2 mod tt, so that 6'i-odd quantities such as ni are discontinuous there. These 
singular behaviors represent the RW phase transition. 



C. Thermodynamics as a function of 6'q and 9i 

Figure [8] presents Q/T^, Im[nq]/T^ and Im[ni]/T^ as a function of 6'q and 6i in the case of 
T = 175 and 250 MeV. The symmetries (|53l) -(l59l) are seen as a bird's eye view. This result 



is consistent with LQCD ones [13711 : in particular, discrete symmetry (1591 ) is clearly seen in the 



LQCD data. If the pion condensate is nonzero, the symmetries (I53l)-(l59l) break down, as shown in 
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SecHni Hence, the fact that LQCD has symmetry (1531) means that the pion condensation does not 
take place also in LQCD simulation. 






Fig. 8: Q/T^, Im[nq]/T'^ and Im[ni]/r'^ as a function of 0q and Oi. Panels (a), (b) and (c) correspond to 
175 MeV, while panels (d), (e) and (f) to 250 MeV. 
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Fig. 9: (a) lm[nq\/T^ and (b) Im[ni\/T^ as a function fo Oq^ and 9i for the case of T = 250MeV. 
D. Comparison of PNJL results with LQCD results 



LQCD data are available at temperatures below and above T^w in Ref- BTII . where the lattice 
size is 16^ x 4 and the forth-rooted KS fermion is taken. The quark and isospin number densities, 
Im[nq]/T^ and Im[?2i]/T'^, shown in Fig. [8] are qualitatively consistent with the corresponding 



LQCD results presented in Figs. 2, 3, 9 and 10 of Ref. BTII . For the case of T > Trw, the 



consistency is more clearly seen in Fig.|9]where Im[r2q]/T^ and Im[ni]/T^ are plotted in the same 



scale, ejTT < 0.3 and Oi/n < 0.3, as Figs. 9 and 10 of Ref. BTfl. In Ref. UM, LQCD data on rzq 
and ni are fitted by the hadron resonance gas (HRG) model 142(1 for the case of T < Tc, since the 
model is one of the most reliable models at T < Tc and also at T = Tc the model is successful 
in fitting the LQCD data by adding correction terms to it. This makes more precise comparison 
possible for T < Tc. 

In the HRG model, Im[nq] and Im[ni] are obtained by sums of free-gas densities over kinds of 
particles [37i]: 



Im[nq]/T=^ = ^ 351^B,7(T)(5(/)sin(3S^q)cos(2M), (60) 

B,I>0 

ImH/T=^ = J2 2/l^B,/(T)5(5)cos(35^q)sin(2/^i), (61) 

B,I>0 

where S{n) = 1 — Snfl/2 and B (I) is the baryon (isospin) number of particle. The parameters are 
fitted to LQCD data in 9q — 9i plane. The resultant values are summarized in Table II. 

Figure [TO] presents lm[nq]/T^ and Im[7ii]/T^ at T = 0.951Tc. The solid (dotted) lines stand 
for the PNJL (HRG) results. In panels (a) and (b) where Im[nq]/T^ is plotted, the PNJL result is 
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W2,l 




O.QSlTc 


0.257 


0.0106 


0.0212 


0.0265 


0.0009 


0.0006 


0.00090 




Tc 


0.3214 


0.0220 


0.0344 


0.0393 


0.0042 


0.0015 


0.0031 


0.010 



TABLE II: Summary of the parameter set of the HRG model in Table IV of Ref. 1371] . 



adjusted to the HRG result at (6'q, 9i) = (tt/G, 0) by multiplying the PNJL result by 4. In panels (c) 
and (d) where Im[?2i]/T^ is drawn, the PNJL resuk is fitted to the HRG result at (6'q, 9i) = (0, tt /5) 
by multiplying the PNJL result by 6. 1 . Oscillatory patterns of the HGM results are well reproduced 
by the PNJL model. 
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Fig. 10: Comparison of the PNJL model with the HRG model for Im[nq]/r^ and Im[ni]/r^ at T = 
0.951Tc = 165 MeV; (a) 6'q dependence of Im[nq]/r3 at 6'i = 0, (b) dependence of lm[nq\/T^ at 
9q = 7r/6, (c) ^q dependence of Im[ni]/r^ at 9i = vr/S and (d) Oj dependence of Im[ni]/T'^ at ^q = 0. 
The soUd (dotted) lines denote the PNJL (HRG) results. The PNJL result is multiplied by 4 to fit the HRG 
resuh at (6'q, 6'i) = (vr/e, 0) in panels (a) and (b) and by 6.1 to fit the HRG result at (6'q, 0i) = (0, vr/S) in 
panels (c) and (d). 



In Fig. [TH the same analysis is made for T = Tc. Again, the PNJL result is adjusted to the 
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HRG result at {9^, 9i) = (vr/G, 0) by multiplying the PNJL result by 2.15 in panels (a) and (b) 
and at (^q, 9i) = (0, tt/S) by multiplying the PNJL result by 3.8 in panels (c) and (d). Oscillatory 
patterns of the HRG results are reasonably reproduced by the PNJL model. Thus, the agreement 
between the two models becomes better in magnitude as T increases. For the oscillatory pattern, 



the agreement 


is reasonably good at both T 


= O.QSlTc and T^. 
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Fig. 11: Comparison of the PNJL model with the HRG model for Im[nq]/T^ and Im[ni]/T^ at T = 175 
MeV (~ Tc); (a) 6^ dependence of lm[nq]/T^ at Oj = 0, (b) Oj dependence of lm[nq]/T^ at Oq = 7r/6, (c) 
9q dependence of Im[ni]/T'^ atOj = tt/B and (d) Oj dependence of Im[ni]/T'^ at = 0. Definition of lines 
is the same as in Fig. [TO] The PNJL result is multipUed by 2. 15 to fit the HRG result at (6'q, 6*1) = (vr/G, 0) 
in panels (a) and (b) and by 3.8 to fit the HRG result at {9q, 0i) = (0, vr/S) in panels (c) and (d). 



The success of the PNJL model for the oscillatory pattern may indicate that the pattern is 
essentially controlled by discrete symmetries of (|9l)-(fT6l). For magnitudes of Im(nq) and Im(ni), 
meanwhile, the PNJL model underestimates LQCD results by a factor of 2 ~ 6. Here we consider 
a possible origin of the discrepancy. In Fig. [I2ta), Im(nq) is plotted as a function of T for the 
case of (^q, ^i) = (7r/6, 0). At T = 1.25Tc = 216 MeV, LQCD data (plus symbol) is larger than 
the Stefan-Boltzmann high-T limit (dot-dashed line). Meanwhile, the PNJL result (solid curve) 
is smaller than the limit at T = L25Tr. The PNJL model is considered to be reliable above T^. 
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Actually, for real quark chemical potential, the PNJL prediction on Uq is consistent with LQCD 
data [17]. We then normalize the LQCD data so that the data at T = 1.25Tc can agree with the 
PNJL result at T = 1.25Tc. The normalized data are shown by cross symbols. At T = 0.951Tc = 
165 MeV and Tc = 173 MeV, the PNJL result is smaller than the normalized data by a factor of 
about 2. This discrepancy is understandable as follows. Below Tc, in general, hadronic excitations 
are important, but such an effect is not included in the mean-field approximation used in the present 
PNJL calculation. The ideal-gas model is considered to be good for T < Tc where hadrons have 
no decay modes. The ideal-gas model yields = 0.0315 for proton and neutron with physical 

masses. Substituting the value for Wxxi2 in (l60l) and adding this correction to the original Im(?7,q) 
of the PNJL model, we have new Im(nq). The new Im(?2q) is plotted by the dashed line up to T^. 
This line agrees with the normalized LQCD data at T = 165 and 173 MeV. 




T[GeV] T[GeV] 



Fig. 12: T dependence of (a) Im(nq) at ((9q, Qi) = (vr/e, 0) and Im(ni) at (6'q, 6*1) = (0, it/5). LQCD data 



are taken from Il37l] . The original values of LQCD data are plotted by plus symbols. The LQCD data are 
normalized so as to reproduce the PNJL result at T = 216 MeV. The normalized LQCD (n-LQCD) data 
are shown by cross symbols. The dashed line is the result of the PNJL density plus the free-gas density; 
as a free particle we take nucleon for Im(nq) and pion for Im(ni). The dot-dashed line represents Im(nq) 
in the Stefan-Boltzmann limit. The LQCD result is multiplied by 0.53 and 0.60 to fit the PNJL result at 
T = 1.25Tc in panels (a) and (b), respectively. 

The same analysis is possible for Im(?2i). Figure [I2tb) presents Im(ni) as a function of T for 
the case of (6'q, 6*1) = (0, 7r/5). At T = 216 MeV, LQCD data (plus symbol) is larger than the 
PNJL result by a factor of 1 .5. fJence the data are normalized so that the data at T = 216 MeV can 
reproduce the corresponding PNJL result. The data thus normalized are shown by cross symbols. 
At T = Tc = 165 and 173 MeV, the PNJL prediction underestimates the normalized LQCD data. 
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Now the pion free-gas density is added to Im(ni), where the pion mass is taken to be 280 MeV 
(the value of the LQCD calculation IsV]). The new Im(ni) is plotted by the dashed line up to Tc. 
The new PNJL result agrees with LQCD data at T = = 165 and 173 MeV. 

As mentioned in Ref. [37], the HRG model works well at T < Tc, but not T > Tc. At T ~ Tc, 
corrections of a few percent to the model prediction are needed. This property is seen also in the 
PNJL result, as shown below. 

Noting that is 6'q-odd, we can find from (ISTI) that 



afc sm{3k9q) . (62) 



fc=0 

The Ofc terms with k > 2 correspond to corrections to the HRG model. Now, we introduce a partial 
sum 

^max 

nqik^i^x) = ^o.k sin(3A;6'q), (63) 

k 

where the ak are evaluated from calculated with the PNJL model. Figure [13] shows 9q depen- 
dence of Uq and nq(A;max) with some values of /cmax^ where the case of = is taken. Panels (a), 
(b) and (c) correspond to the cases of T = 165 MeV (< Tc), 175 MeV (~ Tc)and 185 MeV (> Tc), 
respectively. The PNJL result, riq = nq(/cmax = cxd), is well approximated by nq{k^s_^ = 1) for 
T < Tc, nq(/i;max = 2) for T ~ Tc and Uq^k^g^^ = 10) for T > Tc, as expected. The relative value 
aa/ai at T ~ Tc is 0.11 for the PNJL resuh, while 0.12 for the LQCD result. Thus, the PNJL 
result is consistent with the LQCD result. 



E. Phase diagram in ^i-T plane 

In this subsection, the phase diagram is explored mainly in /ij-T plane, since the phase diagram 



2M. 



in jiq-T plane has already been analyzed for the case of /xi = in Refs. HQ! 

Figure [14] presents T dependence of the absolute value |^ | and the chiral condensate a for two 
cases of /ii = and n/2, where a is normalized by the value ao at T = fiq = fii = 0. Note 
that is an increasing function of T, while the normalized a is a decreasing function of T. 
When 9i = 0, both the chiral and the deconfinement transition are crossover, as represented by 
the solid curves. Their pseudo-critical temperatures are T^^ = 216 MeV for the chiral transition 
and Tc°nf = 173 MeV for the deconfinement transition in the present PNJL calculation, while 
rpx ~ yconf ^ 173 ± 8 MeV in LQCD calculation \M- Thus, the correlation between the two 
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Fig. 13: 6'q dependence of riq and nq(A;max) for (a) T = 165 MeV< Tc, (b) T = 175 MeV~ and (c) 
T = 185 MeV> T^. The case of 6li = is taken. 



transitions is weaker in the present PNJL calculation than in the LQCD simulation. In the previous 
paper pll] . therefore, we introduced the scalar-type eight-quark interaction in the PNJL calculation 
in order to solve this problem; actually, ^ T^"""^ = 173 ± 8 MeV in the PNJL calculation with 
the scalar-type eight-quark interaction. 

For 9i = 7r/2, as denoted by the dashed curves in Fig.[l4l the deconfinement phase transition 
becomes first order, while the chiral condensate hardly depends on T. As shown in (|52l ). the u- 
quark loop contribution to is nearly canceled out by the d-quark one, when 9i = 7r/2. As a 
consequence of this cancellation in H, a has a weak T dependence, while T dependence of <P is 
controlled by the pure gauge part U. The potential U breaks the center symmetry spontaneously, 
when T >Tq = 212 MeV; as shown in Ref. [33], U has two local minima at \<P\ = and ~ 0.45 
for T near Tq, and the local minimum at \<P\ = is deeper than the other only when T < Tq. 
Eventually, T™'^^ nearly agrees with Tq and hence becomes much smaller than in the present 
PNJL calculation with no the eight-quark interaction; i.e., T^"""^ = 212 MeV and = 455 MeV. 
The difference between and T^"^^ is still large, even if the eight-quark interaction is taken into 
account; i.e., T5°"^ = 212 MeV and = 405 MeV. 
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Fig. 14: T dependence of and a normalized by the value ctq at T = /iq = /^ij = 0; as T increases, |^| 
increases but a decreases. The solid (dashed) curves correspond to the case of 6*1 = (7r/2). 

Since two-flavor LQCD data are not available at 6'i = 7r/2, it is not clear whether the large 
difference is realistic. However, it should be noted that LQCD data at 6*: = 7r/2 are available in 
the 8-flavor case [36]. The data show that the chiral and deconfinement transitions are first order 
and ^ j-conf Unfortunately, it is not straightforward to apply the PNJL model to the 8-flavor 
system, since the LQCD data do not present the pion mass and the pion decay constant and hence 
we cannot determine the parameters of the PNJL model. If the PNJL calculation is done without 
changing the parameters from the 2-flavor case, the calculation shows ^ j-conf therefore 
cannot reproduce the LQCD data. This disagreement of the PNJL result with the LQCD data 
is originated in the fact that the correlation between a and <P is weak in the PNJL model. This 
suggests that of the PNJL model should have a direct coupling term such as a^xp*. Thus, the 
nature of the coincidence between the chiral and deconfinement transitions in LQCD, or the origin 
of the direct-coupling term between a and <P in the PNJL model, is an interesting subject as a 
future work. 

Figure [l5] shows the phase diagram of the deconfinement phase transition in Oj-T plane, where 
panels (a), (b) and (c) correspond to three cases of 6'q = 0, 7r/6 and n/3, respectively. The solid 
curves denote the first-order phase transition, while the dashed lines stand for the crossover tran- 
sition. Near 9i = n/2 mod vr, the deconfinement phase transition are first order in all the cases. 
Near 9i = n mod vr, the deconfinement phase transition is first order when 6'q = 0, but crossover 
when ^^q = 7r/6 and 7r/3. The RW phase transition occurs in the area labeled by "RW" between 
the two dot-dashed lines. The dot-dashed line is a boundary of the area and is called "the RW-like 
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transition line " in Ref. Bol . It is a nearly-vertical line starting from point A and is expressed as 
= 7t/2 - 6{T) where 6{T) is defined in Point A is located at (Ta, ^a) = (1.23Tc, 0.4947r) 
in the present 2-flavor analysis, while 8-flavor LQCD data Im show (Ta, 6a) = (1.2Tc, 0.487r). 
Thus, the present result seems to be consistent with the LQCD data. 
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Fig. 15: Phase diagram of the deconfinement phase transition in Oj-T plane. Panels (a), (b) and (c) are 
the cases of 6'q = 0, 7r/6 and vr/S respectively. The first-order (crossover) transition is denoted by the solid 
(dashed) curves. The area labeled by "RW" between the two dot-dashed lines represents the region in which 
the RW phase transition takes place. Point A is located at (Ta, ^a) = (212 MeV, 0.4947r). 

Figure [16] shows the phase diagram of the deconfinement and the RW phase transition in 9q-T 
plane at = 0. The solid lines represent the first-order deconfinement transition, while the dashed 
lines do the crossover deconfinement transition. The dot-dashed lines stand for the RW transition 
line, while point E denotes an endpoint of the RW transition. In panel (a), the present PNJL model 
reproduces LQCD data [6] at finite ^^q. The phase diagram near the RW endpoint (point E) is 
magnified in panel (b). Thus, the RW endpoint is first order in the present PNJL calculation with 
RRW-type U UM; detailed analyses will be made later in Fig. [TSl However, it was second order 
in the previous PNJL calculation BBII with F-type 14 u4\ in which a form inspired by a strong 
coupling QCD was taken for U. Thus, the order of the deconfinement phase transition near the 
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RW endpoint strongly depends on U taken. For comparison, the previous PNJL result is plotted 
together with LQCD data in Fig. [I7l Thus, the present calculation gives better agreement with 
LQCD data than the previous one. In this sense, the present PNJL calculation is more reliable. 
The result of the present PNJL calculation is consistent with a latest LQCD result [|43ll in which 
the order of the RW phase transition at point E is first order for small quark mass, although it is 
second order for heavy quark mass. 




9^/(71/3) 



Fig. 16: Phase diagram of the deconfinement and the RW phase transition in plane at 9i = 0. The 
solid lines stand for the first-order deconfinement transition, while the dashed lines denote the crossover 
deconfinement transition. The RW transition is denoted by the dot-dashed curve. Point E is an endpoint 
of the RW transition. In panel (b), the phase structure near point E is magnified. Lattice data are taken 



from Ref. 



; the pseudocritical temperature at 6'q = is assumed to be 173 MeV determined from LQCD 



calculation of Ref. 



401. 



Finally, the behavior of the RW transition near endpoint E is analyzed more explicitly. Fig- 
ure [13a) shows T dependence of phase if) of the modified Polyakov-loop dlO^^ = n/3 and 
9i = 0. The solid line shows the PNJL prediction with RRW-type U, while the dashed line corre- 
sponds to the result of F-type W. The phase ip is an order parameter of the RW phase transition OBIl . 
Obviously, the RW phase transition at endpoint E is first order for RRW-type U, but second order 
for F-type U. As shown in Fig. [I6](b), there is a meeting point of the solid and dashed lines at 



T = 0.187 MeV, 9^ = 0.93 x n/3 and 9i = 0. This is a critical endpoint of the deconfinement 
phase transition by definition. Figure [TSlb) presents the chiral and Polyakov-loop susceptibilities, 
Xa and x<f, as a function of T at = 0.93 x n/3 and 9i = 0, where RRW-type U is taken. 
The susceptibilities are divergent at the critical endpoint. Hence, the chiral and deconfinement 
transitions are second order at the critical endpoint. 
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Fig. 17: Phase diagram of the deconfinement and the RW phase transition in 9q-T plane at 6i = 0. Here, 
the Polyakov potential of F-type 11411 is taken in the PNJL calculation. See the figure caption of Fig.[T6]for 
definition of lines and lattice data. 
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Fig. 18: (a) T dependence of phase tp of the modified Polyakov-loop at 9q = vr/S and 9i = 0. The 
solid line shows the result of RRW-type U |[l8]. 



while the dashed line corresponds to the result of F-type 



U (b) T dependence of the chiral and Polyakov-loop susceptibilities, Xcr and at 6q = 0.93 x vr/S 



and (9i = 0. 

V. SUMMARY 



We have explored the phase diagram of two-flavor QCD at imaginary quark-number and isospin 
chemical potentials, yUq = iT9q and fiiso = iTO^so- At imaginary /ijso, the pion condensation does 
not take place. The QCD vacuum is then symmetric. As a consequence, at imaginary niso and 
/iq, the partition function (the thermodynamic potential) has discrete symmetries (l9l)-(fTT1) that are 
not present at real /iiso and fi^. The PNJL model possesses all the discrete symmetries, and hence 



the PNJL results are qualitatively consistent with LQCD data presented very lately |136, 



m. In 
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particular, LQCD data BTII have symmetry (|59| ) derived from (l9l)-(fT6l). This indicates that the pion 



condensation does not occur in the LQCD calculation. 



M is made at T < 



A quantitative comparison of the PNJL model with LQCD data pq . 
by using the hadron resonance gas (HRG) model that can reproduce the LQCD data there. As for 
Im[?2q] and Im[rii], the PNJL result underestimates the HRG result in magnitude, but for and 9i 
dependences the agreement between the two is reasonably good. Thus, the PNJL model is useful 
at imaginary /iiso and /iq. 

The PNJL model predicts that the RW phase transition occurs at 6'q = vr/S mod 2n/3 when 
-n/2 - 6{T) < 9i = < n/2 + 6{T), while at ^q = mod 2n/3 when n/2 - 6{T) < 

01 < 37r/2 + 5{T), where 5(T) is given in (|541) . For the case of 9i = 0, the RW phase transition is 
first order at the endpoint in the present PNJL calculation. This is consistent with the latest LQCD 
data l^st]. In a forthcoming paper, we will analyze the relation between imaginary and real 6i. 
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